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Blazars are the most abundant class of known extragalactic very-high-energy (VHE, E>100 GeV) 
gamma-ray sources. However, one of the biggest difficulties in investigating their VHE emission 
resides in their limited number, since less than 60 of them are known by now. 

In this contribution we report on H.E.S.S. observations of the BE Lac object PKS 1440-389. This 
source has been selected as target for H.E.S.S. based on its high-energy gamma-ray properties 
measured by Fermi-LAir. The extrapolation of this bright, hard-spectrum gamma-ray blazar into 
the VHE regime made a detection on a relatively short time scale very likely, despite its uncertain 
redshift. H.E.S.S. observations were carried out with the 4-telescope array from Eebruary to May 
2012 and resulted in a clear detection of the source. Contemporaneous multi-wavelength data are 
used to construct the spectral energy distribution of PKS 1440-389 which can be described by a 
simple one-zone synchrotron-self Compton model. 
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1. Introduction 

Active galactic nuclei (AGN), in particular blazars, are the most abundant extragalactic ob¬ 
jects detected at very high energies (VHE; E > 100 GeV), constituting a third of the known VHE 
gamma-ray sources. Due to their fast flux decrease with increasing energy, observations at VHE 
are mostly performed with ground-based imaging atmospheric Cherenkov telescopes (lACTs). The 
small field of view of current lACTs makes new source discoveries from a large-scale sky survey 
difficult and therefore the search for VHE-emitting blazars has historically involved targeted obser¬ 
vations of source candidates selected based on their radio and X-ray spectra [1]. With the launch 
of Fermi in June 2008, this has partly shifted towards a selection based on extrapolations of high- 
energy (HE; 100 MeV< E < 100 GeV) gamma-ray spectra into the VHE regime [2]. This lead to 
several new blazar discoveries and motivated the choice for targeted observations of PKS 1440-389 
with H.E.S.S. and, eventually, its VHE discovery in April 2012 [3]. 

PKS 1440-389 (R.A. = 14'*43'"57^ deck = -39°08'39", J2000 [4]) is a bright blazar in the 
Fermi-LPH energy regime and has a hard, well-constrained HE spectrum (spectral index: r 2 FGL = 
1.77 ±0.06 [5]). Its distance was estimated as z = 0.065 (6dE Galaxy Survey [6]), but this redshift 
value is not included anymore in the final version of the 6dE catalog due to poor spectral quality 
[7]. Despite many measurements in different wavelength ranges, the redshift of PKS 1440-389 
remains unknown with the currently best constraint to be 0.14 < z < 2.2 [8]. 

2. H.E.S.S. data analysis and results 

H.E.S.S. is an array of five imaging atmospheric Cherenkov telescopes located in the Khomas 
Highland in Namibia which is sensitive to gamma-ray energies from a few tens of GeV to about 
100 TeV [9]. Observations presented here were taken before the installation of the fifth telescope 
for which the array characteristics are given in [10]. 

Observations of PKS 1440-389 were taken between 29 Eeb and 27 May 2012 at a mean zenith 
angle of 17 degrees. All data were taken in wobble mode, for which the source is observed with an 
offset of 0.5° with respect to the center of the instrument’s field of view to allow for simultaneous 
background measurements [11]. After quality selection and dead time correction, the data sum up 
to a total live time of about 12 hours. This data set was analyzed with pre-defined cuts, optimized 
for a steep-spectrum source as decribed in [12]'. The signal was extracted using a reflected region 
background model with an ON region of 0.11° radius centered on the position of PKS 1440-389. In 
this ON region, 183 excess events have been detected, corresponding to a significance of 9.1 a using 
Eq. 17 in [14]. The excess is found to be consistent with a point-like source and the corresponding 
sky map is shown in Eigure 1. A fit to the uncorrelated excess map yields a position for the excess 
consistent with the radio position of PKS 1440-389 [4]. 

The differential energy spectrum of the gamma-ray emission was derived using a forward¬ 
folding maximum likelihood fit [15]. The photon spectrum, shown in Eigure 2, is well-described 
by a power-law function^ with index E = 3.61 ±0.34 and a flux normalization of Nq = (7.98 ± 
1.22) X 10^'^cm^^s^'TeV^' at the decorrelation energy of Eq = 418 GeV (statistical errors only). 

'Results were cross-checked with an independent analysis package [13], yielding compatible results. 
^dN/dE=No-{E/Eo)-'^ 


2 




H.E.S.S. discovery ofPKS 1440-389 


H. Prokoph 



14''5(r 14'’45"’ 


14''40" 

Right Ascension 



14"5<r 14''45'" 14'’4<r 

Right Ascension 


Figure 1: H.E.S.S. sky maps around the position of PKS 1440-389 (marked with a star). (Left) The excess 
map is smoothed with a Gaussian with a width corresponding to the 68% confinement radius of the point 
spread function (PSF) for this analysis (shown inline). (Right) Corresponding H.E.S.S. significance map. 

The integral flux above 220 GeV is (6.81 ± 1.04) x lO^^^cm^^s^^, corresponding to about 3% of 
the Crab Nebula flux [10] above the same energy thi'eshold. A daily binned light curve was derived 
for E > 220 GeV and is shown in Figure 3. A constant fit to these flux points showed no significant 
deviation from a steady flux (x^/ndf = 19.41/14). 
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Figure 2: Differential energy spectram of the VHE gamma-ray emission of PKS 1440-389. Overlaid spec¬ 
tral points were rebinned, requiring a minimum point significance of two sigma per bin. 


3. Multiwavelength data analysis and results 
3.1 Fermi/LAT 

The Large Area Telescope (LAT) on board the Fermi satellite is a pair-conversion gamma- 
ray detector, sensitive in the energy range from 20 MeV to > 300 GeV [16]. Data from Aug 
2008 to Feb 2015 were analyzed using the LAT ScienceTools (version v9r33p0) together with 
the P7REP_SOURCE_V15 instrument response functions. An unbinned likelihood analysis was 
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applied to photon events with energies from 300 MeV to 300 GeV which were selected in a circular 
region of 10° radius centered on the position of PKS 1440-389, resulting in an excess of about 49a 
significance. Two spectral models were applied: a power law and a log-parabola. The results are 
summarized in Table 1. Using a likelihood ratio test, the log-parabola model was preferred with 
respect to the power-law model at the 4.1a level. A three-monthly binned light curve for the full 
time range was computed, resulting in a probability of about 3.4% for a constant fit /ndf = 
39.33/25). The calculation of the excess variance [17] showed only a small hint for variability 
(Fvar = 0.13 ±0.06) in the HE regime. However, no evidence for variability on a monthly time 
scale has been seen within the data set contemporaneous with the H.E.S.S. observations, shown in 
Eigure 3. 


Table 1: Results of the spectral fit to the Eerm/-LAT data. 


Model 

a 

P 

TS 

Elux (0.3-300 GeV) 
[10^^ ph cm^^s^^] 

Eog-parabolic (EP) 

1.32 ±0.04 

0.083 ±0.008 

2403.9 

8.42 ±0.42 

Power law (PE) 

1.79 ±0.03 

— 

2411.8 

9.95 ±0.53 


3.2 Swift/XRT and UVOT 

The X-ray telescope (XRT) on board the Swift satellite is designed to measure X-rays in the 
0.2 — 10 keV energy range [18]. Target of opportunity observations were obtained on 2012 April 29 
(MJD 56046), following the VHE discovery of PKS 1440-389. The analysis of this 8ks exposure 
was performed using standard tools (HEASoft 6.16, Xspec 12.8.2). Data were grouped, requiring 
a minimum of 20 counts per bin, and then fitted with a power-law model. Photo-electric absorption 
was fit using the Wisconsin cross-section [19] with a fixed value for the Galactic column density 
of Atot = 1.08 X 10^' cm^^ [20], resulting in a photon index ofE = 2.64 ±0.05 andAo = (3.80 ± 
0.12) X 10^^ cm^^ s^^ keV^* at 1 keV. It should be noted, that when Atot was left free during the 
fit, a slightly better fit could be obtained with a column density compatible with the total value. 

Simultaneously with the XRT observations, the Swift Ultra Violet and Optical Telescope (UVOT 
[21]) carried out observations in all six filters. Sky-corrected images were taken from the Swift 
archive, and analysis was performed using the uvotmaghist, uvotimsum and uvotsource 
tasks included in the ETOOES software package. Source counts were extracted using a 5" radius 
for V B and U filters and a 12" for the UVWl, UVM2 and UVW2 filters. The background was ex¬ 
tracted from source-free regions in the surroundings. Eilter UVW1 and UVM2 observations having 
four individual exposures, these were stacked prior to aperture photometry as the uvotmaghist 
task showed no variability between the individual exposures. Count rates are then converted to 
fluxes using the standard photometric zero points [22]. The reported fluxes are de-reddened for 
Galactic absorption following the description in [23], with ^(S — V)= 0.103 mag. 

3.3 ATOM 

The Automatic Telescope for Optical Monitoring (ATOM) is a fully automatic 75cm optical 
telescope operated by ESW Heidelberg on the H.E.S.S. site [24]. Operating since 2005, ATOM 
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Figure 3: Light curve of PKS 1440-389 in different wavelength regimes, (top) Nightly (red) and monthly 
(gray) binned fluxes above 220 GeV measured by H.E.S.S. (middle) Monthly binned fluxes measured by 
Fermi-LAH. (bottom) Nightly optical fluxes measured by ATOM in the R band. 


provides optical monitoring of gamma-ray sources and has observed PKS 1440-389 with high ca¬ 
dence in 2012 in the R band. Following the VHE discovery, additional B, V and J band observations 
were obtained. The ATOM data were analysed using the Automatic Data Reduction and Analysis 
Software (ADRAS) which automates standard data reduction and performs differential photometry 
using five known nearby sources as calibrators. A preliminary light curve is presented in Figure 3, 
showing clear variability in the optical regime contemporaneous with the H.E.S.S. observations. 

4. Spectral energy distribution and modeling 

A spectral energy distribution (SED) of PKS 1440-389 was derived using the above described 
multi-wavelength data sets. As no significant variability was detected in the gamma-ray regimes, 
the spectra resulting from the analyses of the full data sets were used. Additionally, archival data 
from the 2MASS [25] and WISE catalogs [26] were used to complete the low-energy part of the 
SED. 

The overall, non-simultaneous SED of PKS 1440-389 is shown in Eigure 4. It exhibits the 
usual two bumps seen in blazars, one at low frequencies, from radio to X-rays, and the other at 
higher frequencies. This SED was modeled with a one-zone leptonic model [27], taking the extra- 
galactic absorption into account [28]. Within the model the low-energy emission is interpreted as 
synchrotron emission from relativistic electrons in a spherical emission region, moving relativis- 
tically along the jet which is closely aligned with the line-of-sight. The high-energy gamma-ray 
emission is produced via Compton up-scattering off the same electron population which produced 
the synchrotron peak (synchrotron-self Compton; SSC). The non-thermal electron distribution in 
the emission region is determined self-consistently as an equilibrium between injection of a power- 
law distribution, radiative cooling, and particle escape. Due to the uncertainty in the distance to 
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PKS 1440-389, the SED modeling was done for two redshifts^: z = 0.065 and z = 0.14. The results 
are shown in Figure 4. It can be seen that the overall SED is moderately well described by the model 
for the used parameter values which are given in Table 2. The parameter values themselves are well 
within the range of those used for previously detected VHE blazars [29] with a relatively low mag¬ 
netic field strength and a large emission region; resulting in a particle dominated jet {Lg/Le < 0.01, 
Lb and Lg being the magnetic and kinetic jet power, respectively). Models with additional exter¬ 
nal radiation fields yield, in general, jet energetics closer to equipartition (Ls/Le ~ 1), but due to 
the lack of additional observational constrains, applying such a more complex model to the SED 
of PKS 1440-389 would result in strongly underconstrained parameters and seems therefore not 
justified. 



Figure 4: Spectral energy distribution of PKS 1440-389 together with the one-zone SSC model for two 
possible redshifts (see text for details). 
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£Be 
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[10^"^ erg/s] 

[103] 

[103] 


[G] 

[10^^ cm] 




[hr] 

0.065 

0.426 

72 

20 

3.2 

0.02 

12 

10 

50 

0.00507 

86 

0.14 

1.038 

78 

20 

3.2 

0.02 

27 

10 

25 

0.01054 

207 


Table 2: SSC model parameters for PKS 1440-389. The columns are the following; z is the assumed 
redshift; Le is the (kinetic) jet luminosity; 71 is the low energy cutoff energy of the electron distribution; 72 
is the high energy cutoff; q is the spectral index of the electron injection spectrum; B is the magnetic field 
strength; R is the emission region radius; E is the bulk Lorentz factor; and Pesc is the escape time parameter 
with fesc = Pesc ' R/c. Additionally, two output parameters are given; Ese is the resulting relative partition 
parameter Ebc = L^lLg, and dfyar is the resulting minimum variability time scale. 


^Redshifts are converted to luminosity distances using a ACDM cosmology with //q = 70kms ^Mpc *, = 0.3 

and fi/v = 0.7. 
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5. Summary and Conclusions 

In this paper, we reported on the H.E.S.S. discovery of the BE Eac object PKS 1440-389 in 
2012. The source was selected as VHE blazar candidate based on its hard high-energy spectrum 
measured by Fermi-LKT, but is also listed in [30] based on its infrared and X-ray properties. The 
average flux in the VHE regime was about 3% of the Crab Nebula flux and showed no significant 
gamma-ray variability. Multi-wavelength data, contemporaneous with the H.E.S.S. observations, 
were used to derive an SED which is well reproduced with a one-zone leptonic model with param¬ 
eters typical for other VHE blazars. 
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